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ABSTRACT 


Simultaneous measurements of waves, and vertical and 
horizontal water particle velocities were made at the 
breaker-line within the surf zone using a capacitance type 
penetrating wave staff, a pressure wave gauge, and an 
electromagnetic current meter. Wave: measurements were also 
made at seaward and shoreward locations. The wave energy- 
MU ty Spectral components were converted to velocity 
spectral components using linear wave theory. These 
computed values compared well qualitatively with the meas- 
ured velocity spectra. Quantitatively, the results showed 
mee linear Theory underpredicted wave-induced horizontal 
velocity spectral components by about 50 percent at the 
frequency of peak energy. The coherence values between 
waves and horizontal velocity were high, ranging above 0.75. 
The phase angle computation showed the calculated velocity 
components leading the measured velocity components by an 
average of 20 degrees, indicating an unstable wave crest 
leading the ie moten ii eae bod or the Wave. 
Probability density functions were computed and compared to 
Gaussian and Gram-Charlier distributions using the chi-square 
goodness-of-fit test. The Gram-Charlier distribution 


gualitatively gave the better fit to the data. 
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A. NATURE OF THE PROBLEM 

The problem is to describe the kinematics of breaking 
waves within the surf zone in terms of the water surface 
elevation. The theories which have been developed for deep 
Water waves can be applied with some: degree of certainty, 
and can be tested in a laboratory situation. These theories 
generally hold up to the time that the wave begins to shoal, 
“rn somewhat lesser accuracy, throughout the shoaling 
Process up to the point of breaking. At the breaker point, 
however, the prediction of wave kinematics becomes more 
mu ftyeult. The type of breaking wave, plunging, surging, or 
Spilling, depends upon deep water characteristics, bottom 
type, and bottom slope, hence requiring a different descrip- 
mun ot the Kinematics within the wave for each breaker 
mc. Physical limitations preclude the construction of a 
realistic laboratory model of the surf zone. And finally, 
the hostile character of the surf zone complicates the 
meine of in situ measurements to the degree that relatively 
little data exist from within the surf zone. 

The most forthright approach to the problem of describing 
the kinematics of breaking waves in the surf zone seems to 
be through direct measurements. Advances in instrument de- 
sign have led to simple, sturdy, measuring devices with 
rapid response time which can sense the small scale as well 


as the large scale features. By making the assumption that 





“ccc Ses are Stationary, this assumption being valid 
at least over a short record of 20 to 30 minutes, the waves 
and wave induced particle velocities ean be dealt with as 
stationary random phenomena, thereby allowing the utiliza- 


tion of standard statistical analysis procedures. 


B. REVIEW OF PREVIOUS WORKS 

The study of the kinematies of breaking waves in the 
surf zone has progressed slowly both due to the difficulty 
in making direct field measurements, and due to the diffi- 
culty in modeling the surf zone in the laboratory. 

Early work by Iversen (1953) involved use of a wave 
channel in which waves were generated using a hinged plane 
oscillating flap. Measurements were made using photographic 
techniques. The length of the channel limited the wave type 
to plunging and surging breakers. Adeyemo (1970) made simi- 
lar laboratory measurements using hydrogen bubble and photo- 
@raphic methods. 

Inman (1956) made field measurements of water particle 
velocity by measuring the drag force under the wave in order 
Benter particle velocities. Miller and Zeigler (1964) used 
Doth acoustic and electromagnetic current meters to determine 
the particle motion in the surf zone. Walker (1969) made 
studies using propeller type flow meters. Thornton (1968) 
and Steer (1972) made use of electromagnetic current meters 


wemmicasure the particle motion. 





CTR OBJECTIVE 

The objective of this research is to study the kinematics 
of water particles in breaking waves within the surf zone. 
Simultaneous measurements were made of the instantaneous 
water surface elevation and two orthogonal water particle 
A e ities at a fixed location in the surf zone. Additional 
wave measurements were made at points seaward and shoreward 
EM he primary location. Estimates of the probability 
Uu ILV functions and power spectra of the wave height and 
En "e velocities are to be made. The applicability of 
linear wave theory is to be measured in computing velocity 
spectral components using the energy-density spectrum of 
DEE aves. The computed velocity spectra are to be compared 


with actually measured velocity spectra. 
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II. STATISTICAL ANALYSIS 


It was assumed that the data to be analyzed represents 
realizations of stationary, random processes. This assump- 
tion allows the use of standard statistical techniques. The 
basic analyses are described in the following sections as 
performed on continuous data of eee jr OU Or 
ease of discussion. The actual analyses were performed 


we ically on records of finite length. 


fee PROBABILITY DENSITY FUNCTION 

The probability density function describes the probability 
that a sample time history record x(t) will assume a value 
between x and (x + Ax), and that this approaches an exact 
meeeaot lity as T, the observation time, approaches infinity, 
and as the defined range, Ax, approaches zero. The proba- 


bility density function is given mathematically as: 


ME 
EE Lim "Es m (a 


Ax>0 T>% T 
Where T, c'r amonanc of time that x(t) falls within the 
range (x, x + Ax). 
The computed probability density function was compared 
to the Gaussian and Gram-Charlier distributions and tested 


for comparison using the chi-square goodness-of-fit test. 
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EE BN density function is given as: 


i ETC 





Pea X) = (2.2) 
am 
where o = standard deviation 
u = mean 
and in standardized form, with z = == ; 
2 
Pa, (2) = T >c pio 
“OT 


The Gaussian distribution is completely described by the 
mean and variance, and has zero skewness and kurtosis equal 
emu nree. 

The Gram-Charlier distribution was shown by Longuet- 
Higgins (1963) to better describe the probability distribu- 
mon Ol the instantaneous water surface elevation, "n, in 
moderately non-linear swell in deep water. MES given by: 


e 2 


| m E 3 
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where 
Id 


Z 


+7 


ee a =) s 


rn LINE 
n n 2 2! 


(Hermite Polynomial) 


skewness of record 


3 
| 


my - kurtosis - 3 


12 





lucen be seen that for zero skewness and kurtosis of three, 
eee curcamoracr Gram—Charlier distribution degenerates to 
the Gaussian distribution. 

To apply the chi-square goodness-of-fit test, the random 
Variable, x(t), was divided into K intervals, called class 
intervals, thereby forming a frequency histogram. The 
observed frequency, fa» is the number of observations within 
cn class interval. The expected frequency, Fis is the 
number of observations which fall in the ith interval of the 
theoretical distribution. The total discrepancy between the 


observed and expected frequency is given by: 
a o (2.4) 


As the value of x* increases from zero, the goodness-of- 


fit decreases. 


BELL SPECTRAL ANALYSIS 
ME Energy=Density Spectrum 
Ei eo cowarianee for a record, x(t), is defined 
mathematically by 
srl 


Ry (1) = EM > x, (t)x, (t+t)dt (2.5) 


| 


The value obtained describes the dependence of the values of 
the data set at one point in time on the values at some 


later time, T, inthe dataset. The function is real-valued and 


E 





a'n Ne mum at vr = 0. The above expression 
implies an infinitely long data set. To deal with the 
Bebo data set of finite length, a lag window is 
ied. The choice of lag window for this analysis was the 
Zen las window, because it has the advantage of no 
negative side lobes, which maintains the coherence function 
between its theoretical values of + 1. It also insures 
hosnumerical resonance problems when computing cross-spectra. 


The Parzen window is given by Bendat and Piersol (1966) as: 


1 - 6G? 4 6 (592 Lom n 
= m m 5 ers bc; x. 
m m 
P(t) = T (240) 
1 43 mE 
2 re | ; T= 5 t lyes’ Th 
0 Vr > T, 


The energy-density spectrum is computed from the Fourier 
transform of the covariance function as modified by the 


Parzen lag window. The spectrum is given by: 


00 


Ze, Seelen "Var (ONUS 


nme f is a particular frequency of interest. 
2. Cross-Spectral Density 
The cross-covariance function describes the general 
dependence of the values of one set of data on the other set, 


in contrast to the auto-covariance function which compares a 
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recora against itself. The cross-covariance function is 
given mathematically as: 
a eer)" 
DES E T 7 x, (t)x,(ttt)dt (OUS 
L 

The quantity R4 5(1) is a real-valued funetion, but it is not 
mese ven function, and does not necessarily have a maximum 
porue at T = O. 

A lag window again is applied to the cross- 
E crlance funetion before computing the cross-spectral 
density. Applying the Parzen window and taking the Fourier 
meanstorm, the cross-spectral density is given by: 


S,o(f) = f POR ¿(me Tar (2.9) 


Since, unlike the energy-density spectrum, the cross-spectrum 
oL an even function, it is expressed in real and 


meecinary parts by: 
S,5(f) = C,5(f) - i Q. (f) | Com) 


The real part, called the co-spectrum, is given by: 


co 


C, Sf) = a P(c)LR (1) + Ri 5(-t) Jeos(2nft) dt (O L) 


The imaginary part, called the quadrature spectrum, is given 


by: Q45(f) = 2 fs P(O)LR,, (1) + RL 5(-1)]sin(2nft)dt (2.12) 
mE 
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eee hase Angle 


Equation (2.10) can be rewritten in complex polar 
Dotation in terms of the moaulus of the cross-spectral 


density and a phase angle, 


Salt) = (S.R (COL 


12! 


where €,5(f) = average angular difference by which the 
cross-correlated components of x(t) lead those of x. (t) in 
each spectral band. The phase angle expressed in terms of 
the co-spectra and quadrature-spectra is given by: 
E 


C12 €f) 


(f) * tan -T (2.11) 


= 


— T < < TI 


H. Coherence 
An expression for the comparison of two random 
variables as a function of frequency is ST DLE 
coherency squared, or simply coherence, and is derived using 
the expressions previously stated for energy-density spectra 
and Dor pumossopeevpPalt density. In mathematical terms, the 
coherence function is given as: 
[Sy 1017 


ATLET < = (2.15) 
51 (1382.11) 


where: 
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DRezesherenee Tunerion has the characteristics that when 
E I n meere completely correlated at 
that frequency; and when E E) - 0, the records are 


statistically independent at that frequency. 
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PRE XPEREMEBNTAL PROCEDURE 


i EXPERIMENTAL SITE 

The experiment was conducted at Del Monte Beach inside 
rey Bay adjacent to Monterey, California, on 12 April, 
1973. The beach profile shown in Figure (1) is one of gentle 
slope. Plunging and spilling breakers are very common at 
Paes site, With spilling breakers being the predominant type. 
Tmhe waves at this location are usually severely directionally 
filtered by refraction due to the geometry of the coastline 


and bay, and approach almost perpendicular to the shore. 


B. INSTRUMENTATION 

The primary measuring instruments were placed on a three 
meter steel tower as shown in Figure (2). The tower was 
guyed so that tower movement and vibration were negligible. 
The instruments were arranged so as to be in the same verti- 
cal plane, with the exception of a small e parallel to 
me Wave crests given to the water current meter to prevent 
interaction with the capacitance probe. 

The water particle velocities were measured with an 
Engineering Physics Company water current meter Model 6130. 
The sea surface elevation was measured with a capacitance 
type penetrating wave staff fabricated locally, and with 
an Interstate Electronics Model DP-200 portable wave recor- 


der with SDP-201 differential pressure sensor. 
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Additional instrumentation was placed seaward and shore- 
ward of the main tower. A pressure wave gauge was installed 
seaward of the main tower by mounting it on a weighted sled. 
The sled was towed into position by attaching a winch cable 
from R. V. Acania. Acania, in her anchored position, heaved 
“nd, pulling the sled into position. The mean water at 
the sled during the experiment was 2.69 meters. 

A Baylor Company Wave Staff System Model: 13528R was 
ieee led 12.9 meters shoreward on a steel tower. This wave 
staff was place in a line with the primary tower and the 
sled so as to be on an orthogonal with the predicted wave 
GELS. 

The capacitance wave staff was designed by L. F. 
McGoldrick (1969) at the University of Chicago as an inex- 
Pemetve, but ultra-linear highly sensitive instrument for 
the detection of capillary waves in the 20 to 30 HZ range. 

It operates on the principal that changing the plate dimen- 
“Sof a Capacitor changes its capacitance Br hence its 
srce on an electrical circuit. The probe was constructed 
of one-eighth inch 0. D. bronze wire; a dialectric of one- 
sixteenth inch wall thickness polystyrene tubing which 

med tightly over the wire; and appropriate electrical 
connections. The sea water in which the probe was immersed 
Gemerstuted the second plate of the capacitor. The probe 

was held stationary on the tower, and as the sea level fluc- 
tuated, the capacitance changed accordingly. These fluctua- 


Ms Were Sensed by a transistorized circuit operated by 
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Poe =contalnea ten volt batteries. The electronics 
section was housed in a watertight case which was mounted 
on the tower. This enabled the connecting leads to be 
relatively short, thereby minimizing wire-to-wire capacitance. 
A static calibration of the system was made by immersing 
Eet in a water tank. The resulting calibration is 
shown in Figure (3). 

Ine Interstate Electronics Pressure Sensor used on the 
Dramary tower and at the sled site is a small, unbonded 
“wn gauge bridge. The sea pressure is coupled by a neo- 
prene diaphragm to a silicone fluid filling the interior. 
me arrangement is such that slow changes, such as tides, 
“ost through the hydraulic filter, and the transducer 
Only senses the rapid pressure fluctuations of waves. The 
high-pass filter response time constant is 48 seconds. 

The Engineering Physics Company electromagnetic current 
Mm? measures two orthogonal components of water particle 
velocity through a range of zero to one and one-half M/SEC 
with a maximum output error of one percent of full scale 
reading. The instrument Jo tuvOnrelecrrmcal time constant 
EunEes, 0.2 and 1.0 seconds. The operation of the unit 
is based on Faraday's principle of magnetic induction. 
Faraday discovered that currents are induced in a conductor 
Ss put in motion through a magnetic field. This 


We t1o0n can be expressed by the vector equation: 
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where E is Che induced electric field, U is Che velocity of 
wecıon, and B is the intensity of the magnetic induction. 

une. sea the induced electric field, E, is linearly pro- 
portional to the amount of fluid flowing through the magnetic 
mera. The magnetic field, B, is most intense near the probe, 
and decreases with the distance from the probe according 

to the inverse square law. Thus, the flow velocity near the 
meme 15 weighted more heavily than that away from the probe, 
with only that water which is within a distance of two to 
n prope radii contributing significantly to the flow 
Maceda voltage. The radius of the probe is 3/8 inch. 

Mne flow meter was recalibrated using the method of Kra- 
pohl (1972) and Steer (1972). This procedure involves oscil- 
lating the probe in a water tank. The equipment arrangement 
is shown in Figure (4). The flow meter probe was mounted 
uNcarprrape whieh travelled back and forth on rails. The 
carriage was driven by a variable speed motor geared to an 
eccentric throw arm. The peak carriage eee was calcu- 
lated to be the tangential velocity of the throw arm. The 


values of maximum carriage velocity, V and maximum measured 


c? 
Velocity, Me were determined for different angular freguen- 
cies. Each pair of electrodes was calibrated by orienting 
them parallel to the flow. 

Figure (5)is a plot of the ratio of measured maximum 
velocity and maximum carriage velocity vs. frequency. The 


instrument's response is shown to be flat out to frequencies 


EN rout 0.5 Hz. AE higher frequencies, the cart became 
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mb addition, waves were generated in the tank 
which rected the Low pattern. This upper frequency limi- 
tation was not deemed a significant factor in measurement 

of Ocean waves in which the significant energy was found to 
be below 0.5 Hz. 

The Model 13529R Baylor Company Wave Staff System is an 
instantaneous water level measuring device consisting of two 
@emewonal half inch, 6x19 IWRC SS wire ropes. Together with 
these lengths of stainless steel wire rope, the transducer 
“ces an electrically linear direct current output that 
is proportional to the amount of wave staff above a short 
circuit produced by the water surface. The working length 
“Staff, through previous modifications, had been 
E ug fron 50 to 5 feet. The results of the recalibra- 


Dare shown in Figure (6). 


C. DATA COLLECTION 
Ehe experiment was performed on 12 April, 1973. The 
towers and sled had been installed at low tide two days pre- 
“Vv. Tests were conducted and final calibration checks 
were Bde cat? April. Included in the calibration 
ne placing of a two volt D. C. reference voltage on 
tne FM tape to be used during data processing. This insured 
that the signal gains could be checked after digitizing. 
During the afternoon of 12 April, swell from a storm off 
the Washington coast began to arrive, SSES a buildup to peak 


MARC coincident with high tide. During the experimer 
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the significant wave height at breaking was about 1 meter. 
The waves were generally of the plunging type, and were 
breaking in the general area of the main tower. The breaker 
point was visually observed to move seaward and shoreward 
Peeording to the incident wave height. 

All data were simultaneously recorded using a Sangamo 
Model 3500, 14 channel FM magnetic tape recorder. Continuous 
M Oring of data was accomplished using strip chart 


recorders and voltmeters. 


D. DATA PRE-PROCESSING 

The data were transcribed using the CS-500 Analog Compu- 
RGO an eight channel Clevit Brush strip chart recorder 
enal verification. It was digitized and recorded onto 
seven track tape using the XDS-9300 computer in conjunction 
with the CS-500 Analog Computer. The digitized records were 
meem@ererred to nine track tape for analysis by the IBM 360 


computer. 
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IV. ANALYSIS OF DATA 


Figures (7) and (8) are typical records showing that 
the breaking waves were roughly triangular in shape and were 
ver peaked. This corresponds to the general observation 
mnm !ungping breakers at or near the breaker point have a 
K Tos vertical front face, and a sloping profile 
toward their trailing edge. The records show small peaks 
between the primary wave peaks. It is thought that the 
small peaks are harmonics of the primary wave frequency. 

The capacitance gauge showed a more sharply peaked wave 
mad the pressure gauges. This was due to the nature of 
the instruments in that the capacitance gauge could sense 


feees rapid changes in the wave slope. 


A. SAMPLING PROCEDURE 

The data were read at a rate of 3 samples per second, 
Dpuesponding to a sample interval of 0.33 seconds. This 
resulted in a Nyquist frequency of 1.5 Hz. The sampling 
pate was considered to be sufficiently high to avoid aliasing 
of energy un oschespoprtilon ot the spectra which was of 
interest Pron observations of the analog record. 

The length of record analyzed was 21 minutes. The maxi- 
mum lag time was taken as five percent of the record giving 
a Spectral bandwidth resolution of 0.00873 Hz. According to 
Blackman and Tukey (1958), this yields 40 degrees of freedom. 


By applying the chi-square distribution, the 80 percent 
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confidence limits for 40 degrees of freedom are found to 
be between 0.77 and 1.38 of the measured power spectral 


estimates. 


meer ROBABILITY DENSITY FUNCTION 
Mae probability density function, pdf, for each wave 
ana Water particle velocity record was estimated using the 
Procedures outlined in Section . For each record, the 
Standard deviation, variance, skewness, and kurtosis were 
computed and are summarized in Table I. Measured distribu- 
tions were compared with the Gaussian and Gram-Charlier 
Hurrbutions using the chi-square goodness-of-fit test. 
Tie various measured distributions were also plotted against 
these analytical distributions. 
Mc [ace Elevation 

At the seaward station, the measured wave pdf shows 
Some relation to the Gaussian distribution [Figure (9)]. 
i aram—Charlier gives a better fit in terms of the chi- 
square goodness-of-fit Lesus ents ss to be expected in 
tnat the Gram-Charlier pdf is calculated using the assigned 
Basmeters from the data, and including moments higher than 
mee second.. 

“mc maey Station, the pdf again follows the 
Gram-Charlier more closely than the Gaussian, but not to as 
great a degree as in the deeper water case [Figures (10) 
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At the shoreward station, the distribution again 
more closely resembles the Gram-Charlier [Figure (12)]. 

The positive skewness measured at each station indi- 
cates a large number of small negative values and a higher 
probability of large positive values of wave height and 
i o9city when compared to the Gaussian. This qualitative 
Seieervation supports the notion that in shallow water, the 
wave crests are steep and narrow, and the wave troughs are 
flat and wide. 

a Water rertiele Velocity 

The probability density functions shown in Figures 
(13) and (14) for water particle velocity again follow the 
Gram-Charlier distribution more closely than they do the 
Gaussian. The Gram-Charlier distribution was computed by 
Carrying the series to include the fourth moments of the 
/ieecuency distribution (Equation 2.3). This appeared to 
be sufficient at the seaward station [Figure (9)]. However, 
BRO the more non-linear situation in the surf zone, it 
appears that higher moments should be included. This would 
Pave a better fit NIE RO hohe mams in the tails of 


the distribution. 


Cee SPECTRAL ANALYSIS 

The energy-density spectrum was calculated for each of 
mm La records. Cross-spectra were then computed between 
Aa ables, from which the coherence and YN were determined. 


Wave energy-density spectral components were computed from 
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pressure energy-density components using linear wave theory. 
Linear wave theory was also used to transform the wave spec- 
I! components to spectral estimates of horizontal and ver- 
tical velocities. These calculéted spectral values were 
Buen compared to measured particle velocity spectra. 

ie Weave Height 

The capacitance wave gauge was taken as the primary 
wave height measuring instrument. The capacitance gauge gave 
a direct measurement of instantaneous sea surface elevation, 
and by taking the mean, a measure of the average sea surface 
elevation was obtained. 

The waves were also measured using a pressure wave 
moe. To spectrally convert these pressure fluctuations to 
mnes of water surface elevation, a transfer function was 
K ed. This function is based on linear theory, and is 


ewen as: 


2 - (El) 
S(f) = ey) S0) 
(wave ) pg cosh k(htz) (pressure) 


H 


where: S(f) energy-density spectra 


k = wave number 
h = water depth 
Z = pressure sensor depth 


The energy-đensity spectra of both the capacitance 
and pressure wave gauges are shown in Figure (15). The 


+. 


pressure spectrum nas been. converted to a Wave speccr uM 
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using Equation 4.1. The spectra, which are characteristic 

of all the records, show a narrow banded energy-density peak 
ae a frequency of 0.07 Hz., corresponding to a period of 15.3 
seconds. Also evident are peaks at 0.14 and 0.21 Hz. which 
appear to be harmonics of the primary peak at 0.07 Hz. These 
harmonics appear to be physical as was observed from the Bu EJ D 
chart earlier. These harmonics probably also have energy 
contributions due to the Fourier computational technique. A 
low frequency peak is also shown at about 0.02 Hz. corres- 
ponding to a period of 50 seconds. Energy at this frequency 
is commonly measured along this coast and is attributed to 
surf beat. 

A peak occurs in all spectra at 0.6 Hz. This is 
attributed to mechanical noise in the tape recorder. This 
peak is at a frequency outside the range of appreciable wave 
ESergy and does not affect the results of the experiment. 

The magnitude of the wave spectral components at 
the frequency of peak energy density calculated from the 
pressure record using linear wave theory was approximately 
60 percent of that measured by the capacitance wave gauge. 
The coherence values were very high, ranging to greater 
then 0.9 in the maximum energy portion of the wave spectrum. 
The phase spectrum shows the capacitance wave gauge measure- 
ments generally leading the pressure measurements by a few 
degrees, probably due to the ability of the capacitance 


gauge to better measure the curling portion of the wave. 


HH 





The spectra of the offshore wave gauge and the capac- 
itance wave gauge were compared, as were the inshore and 
capacitance gauges [Figures (16) and (17)]. A dynamic pressure 
correction [Equation 4.1] was not applied to the seaward 
pressure gauge. The results show that the waves were 
shoaling throughout the process; that is, there was a gradu- 
al buildup of energy as the wave progressed from outside 
the breaker zone into the breaker zone. 

The greater variance at the inshore wave gauge com- 
pared with the wave gauge at the breaker line is attributed 
to energy reflection from the beach. The Baylor wave gauge 
tower was approximately 25 meters from the point at which 
mean water level intersects with the beach, and 12.5 meters 
Shoreward of the mean breaker line where the flow meter and 
primary wave gauges were located. These distances correspond 
to one-half and one-quarter wave lengths, respectively, of 
a 13 second period progressive wave in a depth of 1.9 meters. 
This arrangement places a node at the breaker ee and an 
anti-node at the Baylor wave gauge tower for a standing wave 
Betsch period. This 13 second period corresponds closely 
to the period of peak energy. 

Me phase for a theoretical progressive Wave Mis 
plotted with the measured phase for comparison in Figure (17). 
The phase relationship was calculated using linear wave 
theory. Deviations from the calculated progressive phase 


angle could be attributed to wave reflections. 
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2. Sea Surface Elevation and Water Particle Velocity 


The spectra of sea surface elevation and water 
particle velocities were compared by using linear wave 
theory to convert from wave spectral components to water 
particle velocity spectral components. The water surface 
elevation described by linear theory is given by: 

n = a cos (kx-0t) (4.2) 


where 
sea surtace etevatilon 


a = wave amplitude 

k = wave number 

tea nOorlzonval cartesian eordinare 
o = frequency 


t = time 


The velocity components are given by: 


ao cosh k(h + 2) cos kr come) 


Horizontal velocity = sinh Kh 


l 
C 
H 


sinh k(h +2) _ 
ao Sn Sne a a 


1 
E 
n 


Vertical velocity 


The horizontal velocity spectrum calculated from the wave 


spectrum is given by: 


cosh _k(h +z) 42 


u cos: 4.5) 
S(D (ve) 7 EO “Sinn kh 


ee) 


and the vertical velocity spectrum by: 
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sinh k(h +z) 


: E 2 
E 61) Lo poner SN (4.6) 


where the terms in brackets are the corresponding transfer 
freet on e 
a. Horizontal Water Particle velocity 

Linear theory under-calculates the horizontal 
velocity spectrum by 50 percent at the peoakocenemby mesut 
[See Figure (18)]. The values of coherence are hi 
ranging above 0.75. The phase angle, which according ato 
linear theory should equal zero, ranges up to 43 degrees 
from its average value of about 20 degrees. This shows the 
crest of the unstable breaking wave leading the horizontal 
Memeicle velocities in the body of the flow. This compares 
well with the laboratory findings by Adeyemo (1970) that 
the maximum shoreward horizontal velocity for monochromatic 
waves occurs 25 degrees after the passage of the wave crest. 

Der rtical Water Particle Velocity 

The measured rms vertical velocity was an order 
magnitude less than the measured rms horizontal velocity, 
and showed very little wave induced motion [Figure (19)]. 
The coherence reached a maximum value of 0.37 at a frequency 
EE 5 Hz., which is above the frequency of maximum energy 
measured for the wave and horizontal velocity spectra. 
Meenase angle of 130° to 145° was computed, in contrast to 


the theoretical value of 90° between wave height and vertical 


peMocity. 
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The lack of coherence between water surface 
elevation and vertical velocity could be due to a noisy 
electrical signal from the instrument, a conversion of 


wave-induced motion to turbulent motion, or both. 


De 





The probability density functions for the various records 
compared better with the Gram-Charlier distribution than with 
the Gaussian distribution when tested using the chi-square 
goodness-of-fit test. This result e the observed 
asymmetrical shape of the waves in the breaker zone and 
points out the importance of including higher order moments 
in describing wave phenomena in very shallow water. 

The values of the horizontal velocity energy-density 
spectral components calculated from wave spectra using linear 
We indicate a qualitative, but not a quantitative rela- 
Merienap. In general, linear theory underpredicts the 
Mmeenitude of the velocity spectra. 

The high coherence between waves and horizontal velocity 
shows that most of the motion in the body of the breaking 
Ee is wave-induced. A similar conclusion could not be 
drawn when comparing waves and vertical velocity. 

The spectral difference between water surface elevation 
and horizontal velocity was near that predicted by linear 
theory. The waves lead the horizontal velocity in spectral 
phase on the average by 20 degrees, implying that the 
"curling" crest of the wave arrives prior to maximum water 
particle velocity. 

en summary, the water particle motion in the breaking 


waves seems to be very much wave induced in the body oí p 


D 








u 


yw. Linear theory under-predicts the magnitude of the 










Locity spectra when computing from wave spectra. 
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